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| ABSTRACT 

o 



Context. To investigate the physical nature of the 'nucleated instability' of proto giant planets, the stability of layers in static, radiative gas 
spheres is analysed on the basis of Baker's standard one-zone model. 

Aims. It is shown that stability depends only upon the equations of state, the opacities and the local thermodynamic state in the layer. Stability 
^ . and instability can therefore be expressed in the form of stability equations of state which are universal for a given composition. 
— . o Methods. The stability equations of state are calculated for solar composition and are displayed in the domain -14 < lgp/[gcnT 3 ] < 0, 
' 8.8 < lg e/[ergg~'] < 17.7. These displays may be used to determine the one-zone stability of layers in stellar or planetary structure models by 
directly reading off the value of the stability equations for the thermodynamic state of these layers, specified by state quantities as density p, 
Q ■ temperature T or specific internal energy e. Regions of instability in the (p, e)-plane are described and related to the underlying microphysical 
£j ' processes. 

CZ) . Results. Vibrational instability is found to be a common phenomenon at temperatures lower than the second He ionisation zone. The 
K-mechanism is widespread under 'cool' conditions. 



Key words, giant planet formation - /(-mechanism - stability of gas spheres 

1. Introduction Phenomena similar to the ones described above for giant 

planet formation have been found in hydrodynamical mod- 
In the nucleated instability (also called core instability) hy- els conce rning star formation where protostellar cores ex- 
pothesis of giant planet formation, a critical mass for static plode (Tscharnuter[l987] BaUuch [79881 . whereas earlier stud- 
core envelope protoplanets has been found. Mizuno QSgJJ ies found quas i_ ste ady collapse flows. The similarities in the 
determined the critical mass of the core to be about 12M ffi (micro)physics, i.e., constitutive relations of protostellar cores 
(M e = 5.975 x 10 27 g is the Earth mass), which is independent and pro togiant planets serve as a further motivation for this 
of the outer boundary conditions and therefore independent of study 
the location in the solar nebula. This critical value for the core 
mass corresponds closely to the cores of today's giant planets. 

Although no hydrodynamical study has been available 



2. Baker's standard one-zone model 



many workers conjectured that a collapse or rapid contraction Jn Ms secti(m ±e one . zone model of Baker ( 1966 (> originally 

will ensue after accumulating the critical mass. The main mo- used tQ smdy ±e Csph€id pulsation mec hanism, will be briefly 

tivation for this article is to investigate the stability of the static reviewed The resulting stability criteria will be rewr i tt en in 

envelope at the critical mass. With this aim the local, linear sta- terms of ^ ^ variables> local timesca , es an d constitutive 

bility of static radiative gas spheres is investigated on the basis j ations 

of Baker's ( 1966) standard one-zone model. , i lri££ i . .. . . , c . u - , . 1f 

1 ' Baker ( 1966 1 investigates the stability of thin layers in selt- 

gravitating, spherical gas clouds with the following properties: 
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Fig. 1. Adiabatic exponent Fi. H is plotted as a function of lg internal energy [ergg '] and lg density [g cm 3 ]. 



- energy transport by grey radiation diffusion. 

For the one-zone-model Baker obtains necessary conditions 
for dynamical, secular and vibrational (or pulsational) stabil- 
ity (Eqs. (34a, b, c) in Baker 1966 1. Using Baker's notation: 

M r mass internal to the radius r 

m mass of the zone 

ro unperturbed zone radius 

po unperturbed density in the zone 

To unperturbed temperature in the zone 

L r o unperturbed luminosity 

Zs t h thermal energy of the zone 

and with the definitions of the local cooling time (see Fig.0 

Lrl) ' 

and the local free-fall time 



Table 1. Opacity sources. 
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Baker's K and ctq have the following form: 
n 1 

V32 1 T ff 

7T 6 T co ' 

where « m(fo/po) has been used and 



o- = 



K 



e — mc 



(2) 



(3) 



(4) 



(5) 



is a thermodynamical quantity which is of order 1 and equal to 
1 for nonreacting mixtures of classical perfect gases. The phys- 
ical meaning of cr and K is clearly visible in the equations 
above. <xq represents a frequency of the order one per free-fall 
time. K is proportional to the ratio of the free-fall time and the 
cooling time. Substituting into Baker's criteria, using thermo- 
dynamic identities and definitions of thermodynamic quanti- 
ties, 



Tt = 



V ad = 



(d\nP\ 
\d\npj s 

(dlnT 



v dlnP 



X T = 



d\nP\ 
d\npj T 

d\nP 



d\nT 



din a- 
dlnP 

dh\K 
d\nT 



one obtains, after some pages of algebra, the conditions for sta- 
bility given below: 
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(6) 



(7) 



(8) 



Source 


77 pq 


Yorke 1979, Yorke 1980a 


< 1700" 


Kriigel 1971 


1700 < T < 5000 


Cox & Stewart 1969 


5000 < 



a This is footnote a 

For a physical discussion of the stability criteria see Baker 
fT%6l l or Cox fHWfl . 

We observe that these criteria for dynamical, secular and 
vibrational stability, respectively, can be factorized into 

1. a factor containing local timescales only, 

2. a factor containing only constitutive relations and their 
derivatives. 

The first factors, depending on only timescales, are positive 
by definition. The signs of the left hand sides of the inequal- 
ities and (|8j therefore depend exclusively on the sec- 
ond factors containing the constitutive relations. Since they de- 
pend only on state variables, the stability criteria themselves 
are functions of the thermodynamic state in the local zone. 
The one-zone stability can therefore be determined from a sim- 
ple equation of state, given for example, as a function of den- 
sity and temperature. Once the microphysics, i.e. the thermody- 
namics and opacities (see TableQ, are specified (in practice by 
specifying a chemical composition) the one-zone stability can 
be inferred if the thermodynamic state is specified. The zone - 
or in other words the layer - will be stable or unstable in what- 
ever object it is imbedded as long as it satisfies the one-zone- 
model assumptions. Only the specific growth rates (depending 
upon the time scales) will be different for layers in different 
objects. 

We will now write down the sign (and therefore stability) 
determining parts of the left-hand sides of the inequalities l|6}, 
and (|SJl and thereby obtain stability equations of state. 

The sign determining part of inequality © is 1>Y\ - 4 and it 
reduces to the criterion for dynamical stability 

4 

r, > - ■ (9) 



Stability of the thermodynamical equilibrium demands 

X p >0, c v >0, 
and 

X T >0 

holds for a wide range of physical situations. With 

r 3 _ i = JL*l > o 

pT c v 

n =x P +*r(r 3 - 1) > o 



v ad = 



r,-i 



> o 



(10) 

(11) 

(12) 
(13) 
(14) 



G. Wuchterl and C. Ptolemy: Hydrodynamics of giant planet formation 



3 



Fig. 2. Vibrational stability equation of state S v it>(lg <?, lg p). > 
means vibrational stability. 

we find the sign determining terms in inequalities Q and 
respectively and obtain the following form of the criteria for 
dynamical, secular and vibrational stability, respectively: 

3r 1 -4=:5 dyn > (15) 

1-3/4*, 

-(k t -4) + k p + 1 =: S sec > (16) 

4 

4V ad - (V ad K r + k p ) - — =: S vib > . (17) 
31 i 

The constitutive relations are to be evaluated for the unper- 
turbed thermodynamic state (say (po, 7q)) of the zone. We see 
that the one-zone stability of the layer depends only on the con- 
stitutive relations Fj , V ad , Xj > X P ^ K pi k t- These depend only on 
the unperturbed thermodynamical state of the layer. Therefore 
the above relations define the one-zone-stability equations of 
state Sdyn, S sec and S v it>- See Fig. fj] for a picture of S v n,. 
Regions of secular instability are listed in Table 1 . 



3. Conclusions 

1 . The conditions for the stability of static, radiative layers in 
gas spheres, as described by Baker's ( 1966 1 standard one- 
zone model, can be expressed as stability equations of state. 
These stability equations of state depend only on the local 
thermodynamic state of the layer. 

2. If the constitutive relations - equations of state and 
Rosseland mean opacities - are specified, the stability 
equations of state can be evaluated without specifying prop- 
erties of the layer. 

3. For solar composition gas the /e-mechanism is working in 
the regions of the ice and dust features in the opacities, 
the H2 dissociation and the combined H, first He ionization 
zone, as indicated by vibrational instability. These regions 
of instability are much larger in extent and degree of in- 
stability than the second He ionization zone that drives the 
Cepheid pulsations. 

Acknowledgements. Part of this work was supported by the German 
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ABSTRACT 

Aims. We investigate the thermal evolution of strange stars in the 2-flavor color superconductivity and color-flavor locked phases 
under the influence of deconfinement heating. 

Methods. Due to the spin-down of strange stars, the nuclear matter at the base of the thin crusts dissolves into quarks, releasing 
energy to heating the stars. On the other hand, the neutrino emissivities and specific heat involving pairing quarks are suppressed 
by the large pairing gap in color superconducting phases. Then the thermal evolution equation of strange stars is calculated. 
Results. Deconfinement heating delays the cooling of strange stars considerably. The presence of color superconductivity with 
a large gap enhances this effect. Especially, in the color-flavor locked phase, the stars cannot be very cold at an early age but 
they cool slowly. For the stars with strong magnetic fields, a significant heating period could exist during the first several ten or 
hundred years. In addition, we reckoned that a possible theoretical limit line, which is determined by the competition between 
deconfinement heating and surface photon cooling, may indicate the upper limit temperature that isolated compact stars should 
not exceed. 

Conclusions. Deconfinement heating is important for the thermal evolution of strange stars and is especially determinant for the 
stars in color-flavor locked phase which could show characteristic cooling behavior under this heating effect. 

Key words, stars: neutron-stars: evolution-dense matter-pulsars: general 



1. Introduction 

Cooling simulation based on interior physics is of signifi- 
cant interest for the research of compact stars. In accor- 
dance with nuclear physics, a quark matter core could 
be produced in the interior of compact stars (hybrid 
star) and even strange quark matter stars (strange star, 
SS) may exist. Phcnomcnological and microscopic studies 
have confirmed that quark matter at a sufficiently high 
density, as in compact stars, undergoes a phase transi- 
tion into a color superconducting state, which are typical 
cases of the 2-flavor color superc onductivity (2SC) and 
color -flavor locked (CFL) phases l)Shovkovvll2004t lAlfordl 
l2004h . Theoretical approaches also concur that the super- 
conducting order parameter, which determines the gap 
A in the quark spectrum, lies between 10 and lOOMeV 
for baryon densities existing in the interiors of compact 
stars. Recently, the cooling of hybrid stars with color su- 
perconducting quark cores have been investigated. Stars 
with CFL cores behave s imilarly to ordinary neutron stars 
l|Shovkovv fc Eilisll2002l) . And with a designed 2SC+X 
phase, h ybrid stars could also explain the cooling data 
properly ijGrigorian et aLl boOS - ). In these cases, the ther- 
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mal properties of the quark cores are suppressed by the 
large gap, and the hadronic matter parts of the stars play 
an important role in their cooling history. 

The thermal evolution of SSs has also been extensively 
discussed. In early works, it was generally accepted that 
the surface temperature of SSs should be lower than neu- 
tron stars at the s ame age due to the quark direct Urea 
(QDU) processes dAlcock et all Il988l IPizzocherol Il99li 
|Pagelll992t ISchaab et al.lll99r1) . However, since the elec- 
tron fraction could be small or even vanish, the QDU 
processes may be switched off. The cooling of SSs dom- 
inated by the quark modified Urea (QMU) and quark 
bremsstrahlung (QB) processes c an be slower than neu - 
tron stars with standard cooling l|Schaab et alJll997albl) . 
Of course, a color superconducting phase could occur in 
SSs, a nd its effect on the coo ling of the stars is a significant 
issue. iBlaschke et alJ (|200dl) show that the cooling of the 
stars in the 2SC phase (2SS hereafter) is compatible with 
existing X-ray data but that the stars in CFL phase (CSS 
hereafter) cool down too rapidly, which disagrees with the 
data. However, in those calculations an important factor, 
as described below, is ignored. 

An SS, both in normal phase and in color supercon- 
ducting phase, can sustain a tiny nuclear crust with a 
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maximum density below neutron drip (~ 10 11 g cm -3 ) and 
mass M c < 1O~ 5 M due to the existence of a strong elec- 
tric field on the quark surface ijAlcock et alJll98ri lUsovl 
l2004HZheng k Yuil2006^ . The spin-down of the star makes 
the matter at the bottom of the crust compress. As soon 
as the density exceeds neutron drip, the surplus matter 
in the crust falls into the quark core in the form of neu- 
trons. Consequently, the engulfed neutrons dissolve into 
quarks, and the released energy during th is process leads 
to a so -called deconfinement heating fDEO. lYuan &: Zhand 
( 199§) claim that DH delays the cooling of SSs in normal 
phase (NSS) and may even lead to a slight increase in 
the temperature at the early ages of a star under specific 
conditions. 

We argue that the heating effect on the thermal evolu- 
tion of stars in color superconducting phase is huge com- 
pared to NSSs because of the suppression of specific heat 
and neutrino emission involving pairing quarks. Therefore, 
focus on the effects of DH on the cooling of 2SSs and 
CSSs in this paper. Our paper is arranged as follows. We 
recall neutrino emissivities and specific heat, color super- 
conductivity, and the DH mechanism in Scets.2, 3, and 
4, respectively. The cooling curves and the corresponding 
explanations are presented in Sect. 5. Section 6 contains 
our conclusion and discussions. 



2. Neutrino emissivities and specific heat 

The emissivity associated with the QDU processes d-^uev 
and ue^dv of quarks is l)lwamotolll982|) 



E ( D '~8.8xl0 26 a c 



yl/3 T 6 



9 erg cm sec 



(1) 



where a c is the strong coupling constant, pf, is baryon 
density and po — 0.17 fm -3 the nuclear saturation den- 
sity, and T 9 is the temperature in units of 10 9 K. The 
electron fraction Y e = p e / Pb is small and even vanishes 
at a certain set of param eters of pb, a c and s-quark mass 
m s (|Duncan et al.lll~983j) . It has not yet been considered 
that the contribution to the emissivity by the s-u reac- 
tion, which is suppressed by a n extra factor sin 2 fl r ~ 10 -3 
compared to the d-u reaction l|Duncan et al.lll983^ . where 
6 C is the Cabibbo angle. When the QDU processes are 
switched off due to a small electron fraction (Y e < Y ec — 
(3/7r) 1//2 m 3 a c 3 ^ 2 /64. ), the contribution to the emissiv- 
ities by the QMU processes dq-^uqev and QB processes 
did?.— K]-\Q?. yv dominate. T hese emissivities were also esti- 
mated as l)lwamotolll982() 



;2. 83x10 



19 a 2 



Tg erg cm 



sec 



(2) 



e( B >~2.98xl0 19 ( — ) Tq erg cm" 3 sec 



In order to compute the thermal evolution of the stars, 
we also need the th e specific heat involving quarks and 
electrons written as ljlwamotolll98l iBlaschke etalE oOO'l 



^2.5xl0 20 f ^ 



2/3 



-3 v -l 



T 9 erg cm J K 



c P ~0.6xl0 2 



Y e p b 
Po 



2/3 



T 9 erg cm 3 K 1 . 



(4) 



(5) 



Since the mass of the crust is very small, M c <10 M , 
compared with the total mass of the star, its contribu- 
tion to neutrino emissivity and specific heat can be ne - 
glected llGudmundsson et al.lll983t lLattimer et al.1ll994|) . 
Here we also ignore the neutrino emissivity and specific 
heat due to the photon-gluon excitation, because this ex- 
citation is only important for a temperature higher than 
70MeV iBlaschke et al.ll2000|i . which is much higher than 
the typical temperature in our calculation. 



3. Color superconductivity 

It is widely accepted that the color superconducting 
phase is the real ground state of quantized chromody- 
namics at asymptotically large densities. At a certain 
range of the quark chemical potential the quark-quark in- 
teractio n is attractive, driving the pairing between the 
quarks llAlford et a.lJ ll flflfl ll 99Qt I Afford fc Reddlr2003t 
iRapp et allll998HShovkovvll2004|) . Because of the pairing, 
QDU processes are suppressed by a factor exp(— A/fe^T), 
and QMU & QB processes are suppress ed by a factor 
exp(- 2A/fc B T) for T < T c ~0.4A/fc B l|Blaschke et alJ 
2000). In the 2SC phase, two color states of u and d quarks 
pair, whereas the s quark is unpaired. To be specific, 
we suppose that blue-green and green-blue u-d quarks 
are paired , whereas red u and d quarks (u r ,d r ) remain 
unpaired. As a consequence, the QDU processes on the 
red (unpaired) quarks, as d r — > u r ev 1 as well as QMU, 
d r q r — > u r q r eP, and QB, qi r q2r — * qirQ2rVi / , are not 
blocked, whereas other processes involving paired quarks 
are blocked out by a large pairing gap. Therefore, it can be 
estimated that the neutrino emissi yities in 2SC phase ar e 
reduced by about one magnitude l|Blaschke et all FoOOi. 
On the other hand, the neutrino processes involving all 
flavors are suppressed in the CFL phase by the exponen- 
tial factors. For both phases where the specific heat con- 
tributed by the paired quarks is also changed, we apply 
the formula (Bla schke et aDl2000h 



^sq 



3.2c„ 



/ rp 



2.5 -1.7(|;)+ 3.6 (£) exp(-^) 



(G) 



4. Deconfinement heating 



The effect of DH is determined by the number of neutrons 
engulfed by the quark core, in other words, the variation 
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in the mass of the crust. The total heat released per time 
unit as a function of t is 



ffdec(t) 



1 dM c 
rrib dv 



-v. 



(7) 



release per absorbed neutron, 
m the range <?„^10 — 30MeV 



where q n , the heat 
is expected to be 

l|Haensel fc Zduniklll99i|) . and mb is the mass of baryon 
Assuming the spin-down is induced by the magnetic dipole 
radiation, the evolution of the rotation frequency v is given 

by 



3/c 3 



/x 2 z/ 3 sin 2 



(8) 



where I is the stellar moment of inertia, fi = ^BR 3 is 
the magnetic dipole moment, and 6 is the inclination an- 
gle between magnetic and rotational axes. The mass of 
th e crust M r is calculated by a quadratic function of v 
by iGlendenning fc Weberl lll992h . whose result describes 
the cases with intermedi ate frequencies (< 500Hz) very 
well. IZdunik et ail l|200l[) improve the calculation nearly 
up to the Keplerian frequency using a polynomial includ- 
ing terms of higher o rder in v. The mass of the crust reads 



mg terms ot nigner c 
^Zdunik et alJl200l|) 

M c = M c °(l + 0.24i/f + 0.16f|), 

where v 3 = ^/10 3 Hz, and M° r 
the crust in the static case. 



(9) 

1O~ 5 M0 is the mass of 



5. Cooling curves 

The thermal evolution with DH of a star is determined by 
the equation 



C— — L L 
dt ~ v 7 



H, 



dec; 



(10) 



where C is the total specific heat, L v the neutrino lumi- 
nosity, and L 1 the surface photon luminosity given by 



47Ti? 2 crT s 4 , 



(11) 



where a is the Stefan-Boltzmann constant and T s the sur- 
face temperature. The internal struct ure of SSs can be re- 
garde d as temperature-independent ijGlen fc Sutherland! 
1980), and the surface temperature is related to inter- 
nal temperature by a coefficient determined by the scat- 
tering processes occu rring in the crust. In the work of 
iBlaschke et al.l (UpOO), this relation is gi ven by a simple 
expression T s = 5 x 10" 2 T, but following lYuan fc Zhand 
\ 19991). we apply an accurate fo rmula that is demonstrated 
bvlGudmundsson et al.1 l|l983h. 



T s = 3.08 x lOVOT 



(12) 



where g s ,i4 is the proper surface gravity of the star in 
u nits of 10 14 cm s~ 2 ; or for a recent version see the result 
of IPotekhin et alJ Hi 



In our calculations, to be specific, we consider a model 
of canonical SS of 1.4M at a constant density, which is 




2 3 4 5 6 

log(t/yr) 



Fig. 1. Cooling curves of NSSs, 2SSs, and CSSs without 
DH. 



a very good approx imation for S Ss of mass M < IAMq 
l|Alcock et al.Nl986j) . As used bv IBlaschke et alJ l|2000|) . 
we take Y e = 10~ 5 , a c = 0.25, p = 3p for Y e > Y ec , 
which is a representative set of parameters for which the 
QDU processes contribute to the cooling, whereas Y e = 0, 
a c = 0.15, p = 5po for Y e < Y ec . And we also choose A = 
lOOMeV and q n = 20MeV, the initial temperature T = 
10 9 K, initial period Pq = 0.78ms, and the magnetic tilt 
angle 9 — 45°. The gravitational red-shift is also taken into 
account. Then the effective surface temperature detected 
by a distant observer is T s °° = T s yT— R g /R, where R g 
is the gravitational stellar radius. 

We plot the cooling curves without DH of NSSs (solid 
curves), 2SSs (dotted curves), and CSSs (dashed curves) 
in Fig.l. And the observational data , which are taken from 
Tables 1 and 2 in lPage et al.l l(2004h . are also shown in or- 
der to give the readers a feeling of the position of the illus- 
trative curves in the logarithm — t plane. But we will 
not try to fit the data carefully in this paper. It shows that 
the cooling history of 2SSs are similar to NSSs, whereas 
CSSs become very cold at an early age since the specific 
heat is very small, cooler than 10 4 5 K after lOOOyr for 
Y e = 10~ 5 or several hours for Y e = (this curve is not 
shown in Fig.l). We can see the curves of CSSs are very 
far from the data. Th ese conclusions are also indicated by 
IBlaschke etHl l(2000|) using the relation of T s = 5 x 10" 2 T. 

Figure 2 shows the cooling behaviors of 2SSs with DH 
for various magnetic fie lds (10 8 — 10 12 G). And the analogs 
for NSSs can be seen in lYuan fc Zhand l|l999l) . We can see 
DH delays the stellar cooling considerably. As discussed by 
lYuan fc Zhand l)l999f) , the stronger the magnetic field the 
more rapid the spin-down, and most of the nuclear matter 
in the crust dissolves during an earlier and shorter time. 
For a 2SS with a strong (B — 10 12 G) field and small elec- 
tron fraction (dotted curve a), a distinct heating period 
exists in the first several ten years. And in the cases of 
weak fields (B < 10 10 G), stars could maintain high tern- 
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Fig. 2. Cooling curves of 2SSs with DH for various mag- 
netic fields (curve a: 10 12 G, b: 10 n G, c: 10 10 G, d: 10 9 G, 
c: 10 8 G) and the curves without DH (thick curve). 




2 3 4 5 

log(t/yr) 

Fig. 3. Cooling curves of CSSs with DH for a strong (B = 
10 11 G) and a weak (B = 10 9 G) magnetic field and the 
curve without DH (thick curve). 

peratures even at older ages (> 10 6 yrs). In the following 
paragraph, we discuss in detail how DH induces tempera- 
ture rise and delays cooling. 

We here pay more attention to the situation of CSSs 
because we argue the existence of a marked heating effect 
relative to the reduced emissions. To be clear, Fig. 3 shows 
the cooling curves of CSSs with DH for both a strong 
(B = 10 n G) and a weak (B = 10 9 G) magnetic field. It is 
obvious that the cooling curves are changed dramatically 
by DH. The strong magnetic field induces a rapid spin- 
down of the star at the earliest ages, which could enhance 
the effect of DH to make it greater than the cooling ef- 
fect at the beginning. As a result, the temperature should 
rise due to the surplus heat until the increasing luminosity 



equals the heating effect, L 7 = i/dec, so a net heating pe- 
riod appears at the earliest ages. On the other hand, the 
temperature of the star with a weak field decreases but 
does not rise, due to the relative greater cooling effect at 
the start until the thermal release is compensated for by 
DH entirely: L 7 = i?dcc- In this case, since the confine- 
ment energy deposited in the crust is released slowly to 
heat the star, the star with a weak field can maintain a 
high temperature even at older ages (> 10 6 yrs). To con- 
clude, both the stars with strong and weak fields, after 
several hundred years (the specific value of the time is 
determined by the specific condition of the star), could 
arrive at an equilibrium between the cooling and heat- 
ing effects. From then on, the temperature could only be 
reduced in order to rebuild the equilibrium when Ha ec de- 
ceases with time, so the cooling of the stars is delayed. 
Due to this delay, the curves of CSSs cannot be in conflict 
with observational data as shown in Fig. 3. Since the neu- 
trino emission involving all quarks in CSSs is suppressed, 
the equilibrium discussed above is only determined by DH 
and photon emission, and has hardly anything to do with 
the interior thermal properties of the star. Therefore, the 
dependence of the cooling on the electron fraction is elim- 
inated after the first several hundred years. Going back 
to Fig. 2, we see that the mechanism described above also 
influences the cooling of 2SSs. However, since the equilib- 
rium (L v + L 1 — i?dcc) should involve neutrino luminos- 
ity, which is larger at high temperature whereas smaller 
at low temperature than photon luminosity, the evolution 
of 2SSs may be more complicated than it is for CSSs (see 
Fig. 4 for detail), i.e., the cooling history can be roughly 
divided into a neutrino cooling stage [L v 3> £ 7 ) and a 
photon cooling stage (L v <L 7 ), and the cooling is sen- 
sitive to the electron fraction; and the moment when the 
equilibrium is achieved could be very different for different 
magnetic fields. 

Figure 4 shows the cooling curves of CSSs with dif- 
ferent magnetic fields. We can see the cooling of the stars 
with any field is delayed, just as we find in Fig. 2. However, 
in comparison with Fig. 2, there is a question of why the 
temperature rise of CSSs can be more significant than the 
one of 2SSs. For both CSSs and 2SSs, as discussed in the 
previous paragraph, the reason for the temperature rise is 
that the heating effect is greater than the cooling effect at 
the beginning: Hd cc > L 7 for CSSs and Hd cc > L v + L 1 
for 2SSs. Since the neutrino term is absent for CSSs, the 
initial difference between the heating term and luminos- 
ity of CSSs is much larger than the one of 2SSs with the 
same magnetic field and initial temperature. On the other 
hand, with the rise in temperature, the increase in the 
luminosity is proportional to T 22 for CSSs (see Eqs.(ll, 
12)) but to T 8 for 2SSs with a small electron fraction (see 
Eqs.(2, 3), where the term of photon luminosity is ignored 
since L v ^> L 1 at high temperature). Therefore, 2SSs can 
achieve the equilibrium easily after a comparatively small 
temperature rise, but the magnitude of the rise for CSSs 
needs to be much larger. In addition, the needed magnetic 
field intensity to induce the temperature rise of CSSs could 
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Fig. 4. The same as Fig.2 but for CSSs. The dash-dotted 
curve is the limit line, and the band-like structure is sup- 
posed to reflect the uncertainty of q n . 



Fig. 5. Cooling curves of CSSs with B = 10 9 G for differ- 
ent gaps. The thick curves correspond to NSSs. 



be smaller than 2SSs since the initial cooling effect of CSSs 
is smaller. 

We link the points where the cooling curves turn down 
together as a line (dash-dotted curve) in logarithm T s °° — t 
plane in Fig. 4. The temperature indicated by the line is 
expressed as a power form T s °^ im = 6.4x 10 7 i _1 ^ 4 K. It may 
be well-founded that any other cooling curves of isolated 
stars (to our knowledge), regardless of star models, will 
be below this line due to the high heating effect and low 
cooling effect in CSSs. Hence we reckon that the indicated 
temperature may be the upper limit to what compact stars 
can reach at a given age. We must emphasize that this line 
found from the cooling curves of CSSs is only determined 
by the equilibrium between DH and surface photon cool- 
ing. 

Finally, we present the cooling curves of CSSs for dif- 
ferent gaps in Fig. 5. It can be seen that the cooling curves 
are almost independent of the gap on a very large param- 
eter scale (A~10 - lOOMeV). 

6. Conclusion and discussions 

We have studied the cooling behaviors of rotating SSs in 
the presence of color superconductivity by considering the 
effect of DH. The thermal evolution of SSs is now quite 
different from previous results, because DH can delay the 
cooling, and color superconductivity enhances this effect 
significantly, especially in the CFL phase. For CSSs, the 
previous discussions point out that the specific heat is de- 
termined by electrons since the contribution of quarks has 
been suppressed. This reduction leads to a very rapid cool- 
ing th at disagrees with observational data l|Blaschke et ID 
200o|. However, when we consider the effect of DH, the 
results should imply that the cooling curves could not be 
in serious conflict with the data. We even find it is pos- 



sible that CSSs reach a higher temperature than other 
kinds of compact stars in their cooling history. The limit 

temperature line should illustrate this co njecture. 

To be specific, as pointed out by lYuan fc Zhand 



( 1999J), a temperature rising period could exist at the early 
ages due to the DH with a strong magnetic field. And we 
argue that the presence of color superconducti vity may 
lead to a significant rise. lYnan fc Zha.nel {l999) propose 
that this phenomena may be a signature of the existence 
of SS. In our opinion, if the theory of the color super- 
conductivity is reliable, observing a young and quite hot 
source may be possible, although we also note that an im- 
portant so-c alled brigh t ness c onstraint has been suggested 
recently bv iGrigorianl l)2005|) . who argues that it is un- 
likely that objects with a given age are hotter than those 
already observed. Statistically, this constraint is a good 
finding, but we think that in theory the possibility of the 
existence of young hotter stars still cannot be rejected ab- 
solutely. Of course, for our model, the early evolution also 
may be changed to a certain extent if we consider the for- 
mation of the crust of the star. On the other hand, for 
those stars with weak fields (< 10 10 G), our results show 
that they can maintain a high temperature at older ages 
(> 10 6 yrs). Unfortunately, these older sources with weak 
fields also have been not d etected up to now (see, for ex- 
ample, IPopqv et alJ ((2003) for the list of close-by cooling 
pulsars). To summarize, at the present point of observa- 
tions, there is no evidence of the existence of extra hot 
sources. It may imply that the model needs some further 
improvements. However, it still should be emphasized that 
the various heating mechanisms in compact stars need to 
be given more importance when we talk about the star's 
cooling when using the so-called standard scenario. 
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